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We calculate the decay distributions for inclusive B —> X s ^ and B — > X u lv decays in the endpoint 
region, where radiative corrections are large. The computation is done using effective field theory 
methods. The matching coefficients are computed to 0(a s ), and the anomalous dimensions to next- 
to-leading order. The final expressions for the differential decay spectra include the complete 0(a s ) 
corrections, and sum the leading and next-to-leading Sudakov series. We present results for regions 
of phase space where the shape function can be expanded in local operators, and give the matching 
coefficients of the resulting enhanced non-perturbative effects to order a 3 . We show that moments 
of the shape function are not given by moments of local operators once perturbative effects are 
included, explain why the shape function and its moments satisfy different renormalization group 
equations, and contrast this with the situation for deep inelastic scattering. We show that there are 
large perturbative corrections in the usual definition of the shape function. This renders incorrect 
previous prescriptions for combining radiative corrections with the shape function. 



I. INTRODUCTION 

The inclusive decay rates and differential decay distri- 
butions for B — > Xe t, B — > X s e + e~, and B — > X u lv 
can be computed ^[ 0, 0, in a systematic expansion 
in powers of Aqcd/w& and a s (rrib) using heavy quark 
effective theory (HQET) and the operator product ex- 
pansion. The operator product expansion for the dif- 
ferential decay distribution fails in certain kinematic re- 
gions, but one can still compute a suitably smeared decay 
distribution. In B — > X s j decay, the operator product 
expansion gives an expansion for the photon spectrum 
dr/d_E 7 containing derivatives of <5-functions of the form 
A QCD 6{n) (m b -2£ 7 ), which are singular in the endpoint 
region 2E 1 ~ Mb- Nevertheless, one can still compute 
the differential decay rate in the endpoint region pro- 
vided one smears over a range of photon energies which 
is large compared with Aqcd- Smearing over a range 
of photon energies of order A converts the expansion of 
dr/d-E 7 from an expansion in singular terms to an expan- 
sion in powers of (Aqcd/A)™. 1 In the resonance region 
defined by A ~ Aq CD /Mb, the operator product expan- 
sion breaks down. In this region the invariant mass of 
the final hadronic system satisfies M\ ~ Aq CD , and the 
inclusive B — > X s ^ rate is computed by summing over 
form-factors for a few exclusive modes. 

In the shape function region defined by A ~ Aqcd , all 
the (Aqcd/A)™ terms are equally important. The most 
sine ular terms A^ CD S^ (m 6 -2£L) can be summed into a 
non-perturbative shape function [|| S(E^) that describes 
the photon spectrum in the endpoint region. The sub- 
leading sing ularities, A^ CD 5 ( - n - 1 \m b - 2E 1 ) give rise to 
subleading shape functions suppressed by Aqcd/w& 0, 



1 See, for example, Ref. where smeared moments for B — » Xs-y 
are computed to order o^A). 



and will be neglected here. The final state invariant mass 
in the shape function region is M\ ~ M^Aqcd <S Mg, 
so the final state is jet-like. In this kinematic region, the 
final state is described by a collinear-quark, and the ap- 
propriate effective th eory to use is soft-collinear effective 
theory (SCET) SHIM El 

The shape function region is important for measure- 
ments of B — > X s -f and B — > X u lv decays, because ex- 
perimental cuts needed to eliminate backgrounds restrict 
the measurable region of the allowed phase space. The 
shape functions for B — > X s j and B — > X u tv decays 
are related, so the shape function can be measured in 
one process and then used in the other [TJ ■ This avoids 
introducing model dependence into the analysis of exper- 
imental results. QCD radiative corrections are large in 
the shape function region, due to Sudakov double loga- 
rithms, which need to be summed. The effect of Sudakov 
resummation in B decays, and its effects on the extrac- 
tion of CKM parameters such as V u b has been studied in 
detail previous ly i n a series of papers by Leibovich, Low 
and Rothstein |l3j . 

In this paper, we will study the shape function region 
in B — > X s j and B — > X U £D using SCET. Our results 
are given for arbitrary q 2 , and so can be used for B — > 
X s e + e~ decay. The weak decay Hamiltonian for B — > 
X s e + e~ is considerably more complicated than for the 
other two cases, and is a sum over the different operators. 
The expressions for B — ► X s e + e~ can be obtained using 
the methods of this paper, and are not given explicitly. 

Another region of interest is when the smearing size A 
is large compared with Aqcd, wib>A> Aqcd- In this 
case, one has an expansion in terms of local operators. 
Some of the non-perturbative corrections are enhanced, 
and are of order (Aqcd/A)" 1 , rather than (Aqcd/w;,)" 1 . 
These enhanced non-perturbative corrections will also be 
computed in this paper. 

The SCET renormalization group equations allow one 
to sum the Sudakov double logarithms in the endpoint 
region. The matching conditions are computed to order 
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a s , and the renormalization group evolution to next-to- 
lcading order. The matching conditions and anomalous 
dimensions for the firstpart of our calculation have been 
computed previously 0. We extend the results to in- 
clude all non-vanishing terms for large moments. In the 
regime A 3> Aqcd, the total order a s part of the calcula- 
tion (i.e. not including renormalization group evolution) 
agrees with an existing calculation by de Fazio and Neu- 
bert [l4|- The total order a s contribution is generated 
at different scales. Our calculation shows that the most 
singular order a s terms, of the form ln(l — x)/(l — x) 
are generated at two different scales, and some of them 
are included in the shape function. As a result, shape 
function effects cannot simply be incorporated by convo- 
luting a non-perturbative shape function with the order 
a s decay distributions, as was done in Ref. |14| . and we 
disagree with their results in the shape function region 
A ~ A QCD . 

We study the connection between the shape function 
and its moments. Unlike the case of deep inelastic scat- 
tering, the shape function and its moments satisfy differ- 
ent renormalization group evolution, and are not simply 
related. The differences arise because of the existence 
of the velocity v in heavy quark decays, which couples 
the + and — lightcone components of momentum, and 
because the twist expansion of deep inelastic scattering 
is not valid for heavy hadron decays. This implies that 
moments of the shape function are not related to matrix 
elements of local operators as was previously assumed. 

The outline of the paper is as follows. The kinemat- 
ics and notation are summarized in Section ^ and the 
important mass scales in Section imi The decay distribu- 
tions and hadronic tensors that we compute are discussed 
in Section IIVI The details of the computation are given 
in Section [V] and the relation between the shape func- 
tion and its moments in Section IV1I Some applications 
are discussed in Section IVIIII The reader not interested 
in theoretical details can skip directly to section IVIIII 



II. NOTATION AND KINEMATICS 

The velocity of the decaying B meson is v 11 , and we will 
work in the rest frame v^ L = (1,0,0,0). For B — > X s ~f, 
B — > X s e + e~ , and B — > X u lv decay, the b quark decays 
into a light quark with momentum p and a gauge boson 
with momentum q. The coordinate axes are chosen so 
that the virtual gauge boson is emitted in the ^-direction, 
and the light quark in the negative z-direction. It is 
convenient to define the null vectors n M = (1,0,0,-1) 
and n = (1, 0, 0, 1), which satisfy 



where 



2v t " 



0. 



Any four-vector a 11 can be written as 



ifl+n^ 1 + ia n M + a±, 



a = n ■ a. 



The _L components satisfy n ■ = n- aj_ = v-n- aj_ =0. 
The e tensor is defined so that e 0123 = — e i23 = 1, and 

e'l = e 1 ' v a n p , 

tff = g aP -(nV +n l3 v a ) + n a n 13 . 

The momentum of the b quark in the B rest frame 
is pt — mi,v + k, where k ~ Aqcd is the residual mo- 
mentum, and the coordinate choice is such that = 0. 
Let q 2 be the invariant mass of the virtual gauge boson; 
q 2 = for B — ► X s j. The b quark decays into a light 
quark with momentum p and gauge boson with momen- 
tum q, pi, = p + q. The momentum components of the 
particles are: 



Pb = m b + k+, p b 
q + = m b x, q~ 
p + = mb(l — x) + fc + , p~ 
Pb± = P± = qj 



rrib + k 

rnlx-q 2 
rribx 

o. 



(i) 



For B — > X s j decay, the momentum components are 
given by Eq. for q and p, with q 2 — > 0. 

The kinematics simplifies in shape function region x — > 
1 with to/,(1 — x) ~ k + ~ Aqcd- For B — > X s j, 



P 



rrib(l — x) + k" 



and for B — > Xe + e and B — » X u iv decays, 



= m b , q 
p + = m b (l — p~ 



— q_ 

mt, ' 

2 2 



Instead of decay distributions as a function of x, one 
can also study their moments. For a function f{x) de- 
fined on x G [0, 1], the moments are defined by 



M N (f) = 



dx 



jv- 



'fix). 



The endpoint region x — > 1 corresponds to large moments 
N — > oo with the heuristic relation 1 — x ~ 1/N. For 
smeared observables in the endpoint region with smear- 
ing width A, the relation is N ~ 1/A. The shape 
function region corresponds to taking large moments 
with rrib/N ~ Aqcd- It is also convenient to define 
N = NeJ E = 1.78AT. We will compute all terms which do 
not vanish in the limit N — > oo. The vanishing terms can 
be obtained by computing to higher order in the SCET 
expansion parameter A. 
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m b , n-p 
P = \/p + P~ = \/ijn\ -q 2 )(l-x) 



p + = m b (l - x) 
Aqcd 



N 



TABLE I: The scales fii-m that we will use. fi4 is a scale of 
order Aqcd at which perturbation theory is still valid, fj,4 ~ 
1 GeV. 



III. MASS SCALES 

B decays in the endpoint region have four impor- 
tant scales which are relevant for our computation. The 
largest scale in the calculation is of order the b quark 
mass, and is chosen to be [i\ = m b or /ii = p~ = 
(m 2 — q 2 )/mi 1 . At this scale, p + and Aqcd are infrared 
scales and can be neglected, so the final state quark can 
be treated as massless. The appropriate effective theory 
to use at this scale is SCET, where the quark is described 
by a n-collinear field. The QCD operators match onto 
SCET currents at n±. Our calculation does not simulta- 
neously minimize logarithms of m b /fj, and so the 
choice of [i\ is a matter of taste, /ii = m b has the advan- 
tage that it does not depend on the kinematic variable q 2 . 
We will assume that p~ is of order m b , i.e. that q 2 76 77^. 

The next important scale is the invariant mass of the 
final hadronic system /12 = \fp 2 = p = {m 2 — q 2 )(l — x). 
Below this scale, the invariant mass of the final state 
hadrons is large, and they can be integrated out. At the 
scale fj,2, one integrates out the final hadronic states by 
computing the time-ordered product of two SCET cur- 
rents, integrating out the intermediate collinear quark, 
and matching on to bilocal heavy quark operators. The 



non-locality of the operators is set by /13 



= p+ 



Below 



the next scale ^3, these non-local operators can be re- 
placed by local operators. Finally, large logarithms in 
the matrix elements of the operators can be minimized 
by renormalizing them at the scale ^4 of order the scale 
Aqcd of non-perturbative dynamics. /i4 is chosen to be 
large enough that perturbation theory is still valid. The 
scales are summarized in Table [I] 

There is one important difference from deep inelastic 
scattering, the existence of the scale ^,3. Boost invari- 
ance in deep inelastic scattering forbids the occurrence 
of such a scale. However, in the case of heavy meson 
decays, boost invariance is broken by the choice of rest 
frame of the B meson. Equivalently, the decay ampli- 
tudes can depend on v, which does not exist for deep 
inelastic scattering, and so one can have the scale p + . 

Finally, keep in mind that the ratio of the largest and 
smallest scales /ii//i4 is at most five in B decays. 



IV. THE HADRONIC TENSOR AND THE 
ENDPOINT REGION 

All strong interaction effects for the inclusive B decays 
studied here can be encoded in the hadronic tensor 

W ap = -~Im T a<3 , 

7T 



where 



i, is = —I 



(/) 



B 



T 



(0) 



B 



2m t 



and J^a and are the quark currents mediating the 
b — > uiv and b — > S7 transition, respectively, 



1 



J [u) = uP Rl a b. 

A factor of f /m b has been included in the b — > sj current 
so that it has the same dimension as the b — > u current. 

The hadronic tensor W a/3 can only depend on the mo- 
menta ph and q, or equivalently, on the velocity of the 
heavy quark, w M and the momentum of the light quark 
The most general tensor structure possible is 

W (u, s)a p = _ga0 W M + v <* v W ( u ' s) + ie a ^v p q a W^ u ' s) 

(2) 



+q a q fS W i ^ s) + (q a v + q?v a ) W ; 



(u,s) 



using the convention of Refs. 0> Q for Wi. The coeffi- 
cient of the W3 term is the opposite of Refs. 0, 0] be- 
cause we use the opposite sign convention for eoi23- The 
scalar functions W\ depend on all possible Lorentz in- 
variants that can be formed from the two vectors w M and 



P" 



There are two such invariants, and we will chose 



them to be n ■ p and n ■ p, 

w (f) = W^ f \n-p,n-p). 

The inclusive differential decay rates for the decays 
B — * X s j and B — > X u iv written in terms of the scalar 
functions W\ (n ■ p,n ■ p) are 



dr s 

dx~. 



2 m b T^x^ Uw[ s) - w!f ] - x 7 m b W 5 (s) 



dxe dz dp 2 



l2m b T { u) 

+ (xt.(z - 
+m b (l - 



2(1 - z + p 2 )w[ u) 



xi) 



z + p 2 )(z 



2xAW, 



(3) 



where 



(«) 



.off I 



m J h 



G 2 F \V ub \ 2 m 5 



32tt 4 

5 
1, 







1927T 3 
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and we have denned the dimensionless variables 



2E 1 
mb 



2Et 
m b 



2v ■ p 
mb 



and xg = \ — xg. djr is the coefficient of the b — > sj 
operator O7 in the weak Hamiltonian at the scale fi%. At 
next-to-leading order cy(mb) = —0.311 [Tlj . 

The decay B — > A s £ + ^~ can be treated in a similar 
way, however one has to take into account the presence 
of left and right handed fermions. The expressions of the 
decay rates in this case can be obtained, for example, 
from Ref. [3|. 

The hadronic tensor can be calculated using an oper- 
ator product expansion for the time-ordered product of 
the two currents in Eq. @. This procedure has been 
explained in great detail in f| and the scalar functions 
Wi are known perturbatively to order a s ^4| and non- 
perturbatively to order 1 fm\ 0] . For these results to be 
applicable, the phase space for the decay has to be dom- 
inated by a region where the invariant mass of the final 
hadronic system is far away from zero. Unfortunately, 
for many experimentally accessible observables this con- 
dition is not satisfied. 

In the endpoint region, where p~ ^> p + , or equivalently 
Ep ^S> p 2 , the traditional OPE breaks down. However an 
operator product expansion in terms of bilocal operators 
is still possible. The leading order operator (in the SCET 
expansion parameter A ~ \/p + /p~) is the bilocal opera- 
tor 



0„(* H 



dx~e- lx k+ 



1 

2^ 

b v 5(in ■ D + k + )b v 



b v (0)Y{0,x-)b v (x~) 
(4) 



where Y is an eikonal Wilson line in the n direction from 
x~ to O.The scalar functions W\ {n-p,n-p) can be writ- 
ten as convolutions of perturbatively calculable Wilson 
coefficients with the matrix element of this operator 

W^\n ■p,n-p)=f dk+C { /\n -p,p+- k+,n)f(k+,n), 



of the inclusive decay rates is encoded in the perturbative 
Wilson coefficients Cj . 

As explained above, there are four relevant scales in the 
problem. The region of phase space considered here and 
for which Eq. J5| is valid requires /ii ^> (12 ^ Aqcd- The 
calculation proceeds as follows: (1) Match from QCD to 
SCET currents at (2) Run from /ii to /Z2 using the 
SCET anomalous dimension 71. (3) Integrate out the 
final hadronic states at ^2 by computing the time-ordered 
product of currents, and match onto bilocal operators (4) 
Run the bilocal operators to some common scale fj, ~ 
1 GeV. The matrix element of the bilocal operator is the 
non-perturbative shape function /(fc + ,^) and Eq. JSJ is 
the final result. 

However, if /X3 3> Aqcd, one can proceed further. The 
running in step (4) is then performed down to the scale 
/X3 and then (5) Match onto local operators at the scale 
/i3 (6) Run the local operators down to the scale /X4 and 
compute matrix elements. The matrix elements of these 
local operators are the usual HQET parameters Xi, pi 
etc. which also occur in the OPE for totally inclusive 
processes. 

Steps (l)-(4) will be performed in section[V] while the 
steps (5)-(6) will be performed in section IVll 



(5) 



where 



(B\O v (k+,ri\B) 
2m B 



(6) 



The lightcone distribution function (shape function) of 
the B meson, /(fc + ), encodes all the non-perturbative 
effects of the inclusive decay at leading order in the SCET 
expansion parameter A, and has to be determined from 
experiment. It satisfies 



f(k + ,n) = for k+ < -A 



(7) 



since the matrix element (B\ O v (k + , /1) \B) has no dis- 
continuities for k + < —A. All the remaining information 



V. HA (/) IN THE SHAPE FUNCTION REGION 

P + ~ Aqcd- 

A. Matching from QCD to SCET 

The one-loop matching condition between QCD and 
SCET is given by computing the on-shell matrix element 
in QCD using dimensional regularization to regulate both 
the ultraviolet and infrared divergences, and keeping only 
the finite par t |17| . The required calculations have been 
performed in Q and we will just collect the results given 
there. 

Since both the light quark field as well as the heavy 
quark field are given by two component spinors in the 
effective theory, only three different Dirac structures are 
possible for heavy to light currents in the effective theory 



J 



eff 



Xi 



Tb„ 



with 



r = 1, 



7 



Li 

Tjl 



where Xn is the collinear light quark field. 

In order to facilitate writing left and right handed cur- 
rents, we will often write a fourth Dirac structure, which 
however is related to the previous three 



u 5 

7±7 



ILV _L 



where, as before, e^" = €^ va ^v a np with eoi23 
collinear field is defined as 



-1. The 



Xn{x) 



E 



e-^ x [W^ n ]{x) 



where p is the label momentum which contains compo- 
nents of order 1 and order A. The order A 2 compo- 
nents are associated with the spacetime dependence of 



5 



the fields. W n (x) denotes a Wilson line of collinear glu- 
ons along the path in the n direction from x to oo . This 
Wilson line is required to ensu re g auge invariance of the 
current in the effective theory jlOj . 

Since for the processes we are considering we only need 
left handed light quarks, there are only two currents in 
SCET. Define 



ill = XnPmftv 



eff 



Xn Pr b v . 



The full theory currents mediating the decays B — > 
X s j and B — > X u lv can be matched onto these effective 
theory currents. Using the results of |3| we find 



uPRj^b -> Cx{n-p, {c 2 (n-p, p-i) + C 3 (n ■ p, [i^v^jf , 



sPRia^vb — > Ca(u ■ p, (it) (n^gux-riug^jl X + C 5 (n ■ p, m) (v^n u - i^Tty + ie^ v ) j| 



,-cfr 



(8) 



The Wilson coefficients Ci-s are given by 



Ci(n ■ p,m b , [i) = 1 - as ^^ F J g ( n . P; m&i ^) 4- ln(rt ■ p) 



C 3 (n-p,m b ,^) 



4m 



C 2 {n-p,m b ,n) = I- ® 8 )g(n-p,m b ,^)- ln(n ■ p) 



47T 

a s {fi)C F I 
4tt 1 



1 — n • p 
2 — An- p + {n ■ p) 2 



(1 — n • p) 



^2 



2n • pln(n • p) 



(1 — n • p) (1 — n • p) 2 



n • p 
1 — n • p 



(9) 



C 4 (n -p,m b ,p) = 1 - f ^ F < ff (n • p, m b , fi) + In - 2 Info ■ p) I , 



C 5 (n-p,m b ,[j,) = 1 - as ^ Cp J g ( n .p iTOb)jU ) -f i n iL + ln(n-p) ( 4 " ^ — ^ 

47r TO b \ 1 — n • p 



where 



p(n-p, m b ,/i) 



1 In 2 + | In - 2 In hl(n ■ p) + 21n 2 (fi ■ p) + 2Li 2 (l-n • p) + + 

2 mr 2 mr mi 12 



6 , 



and 



B. Renormalization of the SCET Currents 



- ~ n ■ p q" 
n ■ p = = 1 -. 



m b 



The matching coefficients for a^ v have an extra ln/i 2 /ml 
contribution which corresponds to the anomalous dimen- 
sion for this operator in full QCD. 

The matching scale can be chosen to be either fi± — n-p 
or /ii = m b At these scales the matching coefficients are 
given by Eq. © with 



g(n ■ p,m b ,n ■ p) = 5 ln(n • p) + 2Li 2 (l-n • p) + — + 6, 
g(fi ■ p,m b ,m b ) = 2 ln 2 (rl • p) + 2Li 2 (l-n • p) + — + 6. 



The running of the currents in the effective theory 
have also been calculated in @- The currents are multi- 
plicatively renormalized, and satisfy the renormalization 
group equation 

M^Cj(^,ra-p) = 7O, ri • p)C l (p,fi ■ p) . 

The leading order (LO) and next to leading order (NLO) 
anomalous dimensions are 



7lo = 



7NLO 



q s Q)C F ^ / 
a a [p)C F 



2tt 



Mi 

| + 21n 
2 \n ■ p 



Mi 



-2B°ZW F In ( — 
(2tt)2 Vmi 



(10) 



6 



where llE 



B 



r l 67 n 
Ca 1 18 _ T 



—nr. 
9 1 



(11) 



Note that the value of /ii in these anomalous dimensions 
is chosen to coincide with the choice made for the match- 
ing scale in Eq. |JSJ. The [i\ dependence in the a s term 
cancels between 7lo and 7nlo- The difference between 
the two possible choices for \i\ in the a 2 term is NNLO. 

The solution to this renormalization group equation is 
given by 



Ci(n) = Ci(^i)exp 



with 



r o( z ) 

n(z) 



4nC F 



where 



PI 




CfPi \, 


PI 




c F 


r5 


Po 




2C F 


B r 


PI 







r o( z ) 



1 +lnz 



Hi 



■n(z) 



J 1 ' v 
1 — In z 



In ; 



(12) 



a. s {m) ' 



r U r 2 
Po = -j Ca ~ 3"^' 



34. 

y 



01 = ~—C\-^C A n s 



10 

y 



2C 



F n f , 



d/x 



Po 2 Pi 3 



Ca = 3, C F =4/3 and n/ = 4 is the number of light 
quark flavors. 



C. Matching to bilocal Operators 

Below the scale fi 2 — p the final hadronic state is heavy 
and can be integrated out. This is done by matching the 
time-ordered product of two SCET currents onto bilo- 
cal operators by integrating out the n-collinear fields at 
the scale fi2- The tree level time-ordered product graph 
is shown in Fig. ^ The tree level graph for the bilocal 
operator is shown in Fig. [5] and matches on to a bilo- 
cal heavy quark operator, the shape function operator 
defined in Eq. The Feynman rules for matrix ele- 
ments of this operator are given by taking the discon- 
tinuity of the diagram, since the operator is a product, 
not a time-ordered product. Spinors in SCET have only 
two components, which can be used to simplify the re- 
sults significantly. One can easily show that after taking 





FIG. 1: Tree graph for the product of two currents. The 
dashed line depicts a collinear particle. 




FIG. 2: Tree level matrix element of the shape function op- 
erator. The double line between the two currents depicts a 
Wilson line in the n direction. 



matrix elements between pseudoscalar B mesons that the 
only non- vanishing time-ordered products have two iden- 
tical currents jf . Furthermore, perturbative corrections 
in SCET are independent of the Dirac structure of the 
currents. As a result, one can show that the only non- 
vanishing time-ordered products in SCET can be written 
in the form 



B 



T 



eff/ 



B 



<£xe~ lq - x ( B 



1 



B 



= -(B\T(p+)\B). 

Combining this result with Eqs. J3J and JHJ we find for 

the decay B — > X u lv 



(u) 



C*(B\W(p+)\B) 



W. 



(u) 



2m B 

,,2 ' C l ' 



n ■ P ~ 1 ^,2 

1 T I T " 

n ■ p \n ■ p 2 



C 2 



c. 



(B\W(p+)\B) 
2m b 



W. 



(«) 



1 



('«) 



where 



2rribn ■ p 
1 

m 2 (n ■ p) 2 
1 



2 (B\W(p + )\B) 



2m B 



2mb(n ■ p) 2 
-C^C^n ■ p 

W(p + 



{cl _ cl) (BW)\B) 
2m B 

C\(2-n-p)-2Cl 
(B\W(p+)\B) 



2m,} 



(13) 



• Im T(p + ). 
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For the decay B 



wk 



wk 



(» 



-> W,7 we find (using q 2 
Ct(B\W( P +)\B) 



0) 



4 
0, 



2m B 



Cf (B|W(p+)|B> 



2771;, 
1 

Am 2 



2m B 



<B|W(p+)|B> 



2777 B 



C| (B\W(p+)\B) 



2mj 



2m B 



(14) 



Gauge invariance implies (using q 2 — 0) that W : 



(») _ 







and 2VF-P" 1 = TObW 7 ^ 5,1 , so these relations are true to all 
orders in perturbation theory. For the photon spectrum 
in B — > X s j, we will need the linear combination 



= 2m b {AW[ s) - W, 



obtained from Eq. © with 

(s) 

this linear combination of W„ is 



-m 6 W 5 (s) ), (15) 
1. The matching for 



W\ s) = m b C. 



(B\W(p+)\B) 
2m B 



(16) 



Comparing the tree-level matrix element of the time- 
ordered product shown in Fig. ^ with the tree-level ma- 
trix element of the shape function operator shown in 
Fig. 121 one finds 

W(p + ) = O v (p+)+0(a s ). 

The one-loop matching condition is evaluated by com- 
puting the graphs in Fig. and subtracting from it the 
one-loop matrix element of the tree-level operator W(p + ) 
in the effective theory, shown in Fig. 0] The matrix ele- 
ment of O v (p + ) at one loop exactly reproduces the graphs 
in Fig. EUa-d), since the coupling of soft gluons to the 
n-collinear quark is identical to the coupling to the Wil- 
son line Y. The matching condition is therefore given 
by graphs Fig. EHe-g). The Feynman rules for graphs 
Fig. EUc-g) do not depend on the external quark field, so 
the graphs have the same value as the matching graphs 
for deep inelastic scattering [l^, with the replacement 
p — > 0, q — > p. The discontinuity of the graph extends 
over the infinite range < p + < oo, rather than a finite 
interval such as [0,1]. For this reason, it is convenient 
to regulate the singularity at p + — by a modified +- 
distribution, called the /j-distribution since it depends on 
the scale /i. The /^-distribution, its properties, and rela- 
tion to the +- distribution are given in Appendix 1X1 The 
resulting expression for W(p + ) is 

W(p + ) = J dr+C(p+ -r+,n)O v (r+), 



C<°>(g+ aO = <%+), 

"Mg+AO %+) 



c (1) 0z + ,m) 



41n^-3 



2 In 2 — - - 3 In — - + 7 - tt 



5(q + ). 



(17) 



D. Renormalization of the Shape Function 



We can study the renormalization of the shape function 
by computing the on-shell matrix element of the shape 
function operator O v (r + ) in a heavy quark target, with 
residual momentum k = 0. 

The tree level matrix element is given by the graph 
shown in Fig. The spin-averaged matrix element is 



/tree = Disc- 



1 



2tt H 
S (r+) . 



7.0+ 



(18) 



The loop corrections to the shape function matrix ele- 
ment are the graphs in Fig. 01 and wavefunction graphs. 
Graph Fig. 0Ja) gives 



In 



l 2 C f 



x Disc- 



csc(2e7r)r(e)/j 2<: 
' (-r+ - 



(r+ + ?,0+ 
In terms of the /^-distribution, 

9 2 C F 



-2c 



la = - 



16tt 2 



1 



~8{r+)~ 






+ 4 











12 ) « r+ > 



ln(r+//t) 6(r+) 



(19) 



The spin-averaged loop graph Fig. EJb) is 
1 



a 2 C 

i^-r 4Diw | - - 2 In (- r - - /'<)' 



16tt 3 

F 



i0^ 



9 2 C 



■<7T- 



-5(r+)-2 

€ 



(r+) 



(20) 

The 1/e is an infrared divergence. 

The heavy quark wavefunction graph vanishes on-shell 
when evaluated in pure dimensional regularization, be- 
cause the integral is scaleless. One can verify that the 
graph has the structure = 1/euv — 1/ £ IR with exactly 
the right coefficient to convert the infrared divergence in 
Eq. I|20|) into an ultraviolet divergence. 
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operator: 



(a) (b) 
FIG. 4: One loop correction to the shape function operator. 



The wavefunction renormalization of the internal Wil- 
son line is zero since n 2 = 0. 

The one-loop matrix element of the shape function op- 
erator is given by the sum of the tree graph, twice graph 
(a), graph (b) and the wavefunction graphs. From the 
infinite parts, we see that the renormalized operator O v 
is related to the bare operator by a convolution with the 
renormalization coefficient 



Of) (r+) 



+ \ - 



di + Z(r+,£+)O v {t 



(21) 



where is the bare operator, O v (r + ) is the renor- 

malized operator, and 



Z(r+,£+) = S(r+-£+) + 




The crucial difference from the renormalization of the 
deep inelastic structure function is that Eq. (|2*2)l does not 
have the restriction l + > 0, so the ^-function is satisfied 
for — oo < l + < t + , rather than < £ + < r + . 

Differentiating Eq. (|21[1 with respect to /i gives the 
renormalization group equation for the shape function 



dl* 1 (T*,l+)O v (l*),{23) 



where 



7 (r 



a s Cf 



(24) 



using Eq. (|22|) . The shape function f(k + , /i) satisfies the 
same equation, since it is the matrix element of O v (k + ). 

The convolution in the renormalization group equation 
Eq. (|23|l makes this equation difficult to solve and so far 
no solution exists in the literature. Methods that are 
used to solve a similar equation for the parton distribu- 
tion functions do not work in this case because of the 
difference in the ^-function mentioned above. 



E. Expressions for 

In this section we combine the results of the previous 
sections and give the final expressions for the scalar func- 
tions Wl , from which all differential decay rates can be 
obtained. As before, we define 

AJ) 



W t (n ■ p,n- p = p + ) = 
J 6k + c\ f \fi-p;p^ 



k+,fi)f(k+,fi), (25) 



where the shape function /(fc + ,^) is given by Eq. JSJ. 
The expressions for (n ■ p,p + — k + , fi) at fi = /12 to 
first order in a s can be obtained from Eqs. ijT3|) . ijTHl . 
©. As explained earlier, one can chose the matching 
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scale from QCD onto SCET to be any scale fj,± ~ m b 
and in Section fV Al we gave results for both p,i = raj and 
Hi = h ■ p. In this sections we give all results with the 
choice ji\ = m b . Wc find 



S(ji2,m b ) 



U)nt \ , oi s (m h )C F (/) 
i\ C{w, fj, 2 ) H " 



4tt 



by>5(w) 



(26) 
O(al), 



where the coefficients a\ , 6 2 are functions of n -p and 
to;,. The function C(w,fj,2) is universal and does not de- 
pend on the decay process. It is given in Eq. (fTTjl . 

The scale factor S(n2,m b ) is the running of two SCET 
currents from the matching scale \x\ — m b to the scale 
/.t = /12 where the OPE is performed. We can write 



S{n2,m b ) = exp2 



rp(z) 
a s (mb) 



where ro(z) and r\{z) are given in Eq. (|12|) (one 
should pick /ii = m b in these expressions), and z = 
cts{l^2) / oi s {m b ) . The remaining coefficients are different 
for the decay B — > X u lv and B — ► X s 7. For _B — * ATufi' 
we find 



> 



> 



1 

4' 



n ■ p 



1 



2m b n ■ p 
0, 

1 



2m b 7l • p ' 



where 



n ■ p 



The coefficients 6 l - lt ' ) are given by 



n • p 

TOf, 



(" 



d " 



-7>9{n-p,m b ,m h ) - \ogn-p — 

I 2(1 — n ■ p) 

2g(n ■ p,m bl m b ) logn ■ p 



n ■ p 

g(n ■ p,m b ,m b ) _ 
m b n ■ p 
2 

m 2 h ■ p(l — n ■ p) 
g{n ■ p,m b ,m b ) 



6- _ 
n ■ p 

log n ■ p in ■ p — 2 
m b n ■ p 1 — n ■ p 
log h ■ p 1 — 2n • p 



m 2 n ■ p (1 — n • p) 2 



1 



Tfi b Ti ■ p m b n ■ p(l — n ■ p) 

log n ■ p (in ■ p — 2) (n ■ p — 2) 



m b n ■ p 



(1 — n ■ p) 2 



with 

g(n ■ p,m b ,m b ) 



21n 2 (fi-p) +2Li 2 (l-n-j3) 



7T~ 
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For the decay £? — > n ■ p — m b . This simplifies 

the resulting expression significantly and we find 



.0 



1 

4' 
0, 

1 



2TO6 ' 

3 



4m? ' 
1 



For the coefficients b 



6« = -3- 

4 S) = o, 



7 



2TOft ' 
TOf,. 



(27) 



24' 



7T" 



6l s) = — 



10 



12m b ' 

^2 



rn 



6 W = _ 

7 



6 
12 



12toi, ' 

r2 ^ 



7T 



TOf,. 



(28) 



The above equations give C^\n ■ p,w,^2), and so 
must be integrated with the shape function /(fc + , /i 2 ) also 
renormalized at the scale H2- One can instead convert to 
the shape function at some other scale (such as ^3) by 
using the renormalization group equation Eq. I|24|) for the 
shape function. 



VI. 



W} f) IN THE PHASE SPACE REGION 

p + > Aqcd. 



As explained earlier, in regions of phase space where 
p + Aqcd an additional matching can be performed at 
the scale /13 = p + . Since the non-locality of the bilocal 
operator O v (r + ) was set by the scale p + , the operators 
below the scale /13 will be local operators. In this section 
we will match onto these operators at one loop and then 
determine their running. 



A. Matching onto local operators 

The bilocal operator O v (p + ) below ^ can be expanded 
in terms of local operators O m , 



O v ( P + ) =J2v m (p+, f x)O r , 



(29) 
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with coefficients V m (p + , /i) , where (at tree ievei) 

O m = b v (in-D) m b v . 

The matrix eiements of O m are of order Aq Cd . The 
coefficients T> m of are of order (p + )~ m ~ 1 . The subleading 
terms in the endpoint region have coefficients suppressed 
by powers of p + jm^. 

The focal operators O m are related to the moments 
of O v (p + ). The matrix elements of the first few local 
operators are given by 

(B\b v b v \B) = 2m B , 
{B\b v {in ■ D)b v \B) = 0, 

2mg 



(B\b v (in ■ D) 2 b v \B) 
(B\b v (m ■ Dfb v \B) 



3 

2mg 



Ai(m) , 



to leading order in l/m^. Xi and p\ are given by the 
matrix elements of the kinetic and Darwin terms 

(B\b v D 2 b v \B) = -AmsX^p), 
gv fj {B\b v (D a G al3 )b v \B) = -4m B pi(/x). 

Since the scale /i3 is now a large scale, the matrix el- 
ement of the operator O v (r + ) can be calculated pertur- 
batively, using on-shell partonic states with residual mo- 
mentum k = 0. At tree level the matching Eq. I|29|) is 
trivial since O v (p + ) = b v 5(in ■ D + p + )b v + 0(a s ), and 
we find 



V m {p A 



-(— 

ml \dp + 



6(p + ) + 0(a s ) . 



The order a s matrix element is given by the finite parts 
of Eq. (EEJ), twice Eq. dJ) and Eg. ffity . 



(b v \O v (p+)\b v )=6(p+) 



a s (p,)Cf 

47T 



5(p+) 



0(p+) 



p- 



ln(p+///) 9{p+) 



(30) 



Since k + — 0, the matrix elements of O m vanish except 
for m = 0, so we get 



V (p + ,fi) = S(p + ) 



+4 



a s C F 

47T 



S(p+) 



\6{ P +y 


+ 8 


rin( P +/ M ) 6{p+y 




[ P + 




[ P + 


i' 



(31) 



This matrix element gives the parton level value for the 
shape function, 



/part(^ + : A 4 ) 



7T 

"g 



6(k + ) -4 



6»(/c+) 



ln(fc + /At) 6»(fc+) 
fc+ 



(32) 



To determine the coefficients T> m for m > requires 
evaluating the matrix clement Eq. 130fl for non-zero k + , 
with 2fc • v = k + + k~ = 0, so that the quark is still 
on-shell. It can be determined, following the arguments 
in Ref 0, in terms of the lowest order coefficient. The 
matrix element Eq. I|30|l for non-zero k + is given by 
the replacement p + — > p + + k + , which follows from 
reparametrization invariance [20j , so that 

(b v (k + )\O v (p+) \b v (k+))=V (p+ + k+, f i) 
(k+) m 
rn 



The matrix elements of b v (in ■ D) m b v are their tree level 
value (k + ) n , since loop graphs are scaleless and vanish, 
so 



Turnip ,M) = 



1 



d 
dp 



V (p + ,ri. 



(33) 



The matching condition Eq. (|3*3^l is evaluated at a scale 
/i ~ p + to minimize logarithms. At one loop, one can also 
generate four-quark operators; the matching coefficients 
for these operators is not computed here. 



1. Local Operators in Deep Inelastic Scattering 

In deep inelastic scattering, the on-shell matrix ele- 
ment of the parton distribution operator vanishes, so the 
analog of Eq. (|32|l is f' 1 ) — 0. An on-shell external state 
has k + ^ 0, k" = 0. The momentum p + in the operator 
and the momentum k + of the external state enter the 
loop integral over the + component of the loop integral, 
but the — component of the loop integral is scaleless, and 
vanishes. In the B meson shape function matrix clement 
Eq. H30|) . the heavy quark propagator v ■ i mixes the + 
and — components of the loop momentum £, so the scale 
p + enters both loop integrals, and the l~ integral is no 
longer scaleless. 



B. Running of the local operators 

Since the matrix elements of the local b v (in- D) m b v op- 
erators are equal to the usual HQET operator matrix ele- 
ments that define Ai, pi, etc., the renormalization group 
equations for the matrix elements of these two sets of 
operators must be identical. One can check that the run- 
ning of the local b v (in ■ D) m b v operators reproduces the 
known renormalization group evolution for Ai and p\ pf| 
for m — 2,3, respectively. There is an important differ- 
ence here between the evolution of the local b v (in- D) m b v 



11 




(a) (b) 
FIG. 5: One loop corrections to the local operator. 



a 

3 
3 
3 

3- 



3 
3 
3 
3 

.3- 



FIG. 6: One loop running of the local operator. 



operators and the corresponding evolution of local opera- 
tors in deep inelastic scattering. In deep inelastic scatter- 
ing, the evolution of the bilocal parton distribution oper- 
ator is equivalent to the evolution of the local twist-two 
operators; the twist- two anomalous dimensions are given 
by the moments of the Altarelli-Parisi evolution kernel 
for the parton distribution. This is not the case for the 
shape function operator and its moments. The shape 
function anomalous dimension at one-loop depends only 
on C'f, whereas the anomalous dimension of p%, which 
is the same as the anomalous dimension of the Darwin 
term, depends on Ca PH - 

The two diagrams shown in Fig. contribute to the 
renormalization of the local operators O m . For the first 
diagram we find for the forward matrix element (k' = k) 



In 



-ig 2 C F 



d d £ [n-{k + £)} r 



9 2 C F 



(27T) 



£} 2 P 



(n • k) m , 
(34) 



where we have only kept the divergent part of the graph 
and used the equation of motion k ■ v = for the external 
states since we are only working to lowest order in l/m&. 
The graph Fig. |5f b) gives 



I b = ig 2 C F ^ 



■E 

i=i 



d d £ r [n.(k + e)] m - j 

(n ■ k) Yin = 0, 



(27T) 



v ■£ £ 2 



since all the divergent terms are proportional to powers 
of [k ■ v), and vanish. The wavefunction graph gives 



have twist one, and two quark operators with twist two 
cannot mix with four-quark operators which have twist 
four, so the box graph Fig. |H| does not contribute to the 
anomalous dimension. In B decay, there is no analo- 
gous twist expansion, and the mixing graphs of Fig. © 
are present. 2 They are precisely the graphs computed 
in Ref. |2l| which give the anomalous dimension of the 
Darwin term whose matrix element is p\, and allow the 
dimension six operator b v (in ■ D) 3 b v to mix into the di- 
mension six operator b v b v qq. 



C. Expressions for W i 



( f) 

The final expressions for W( in the region A ^> Aqcd 
are given by combining Eq. i|25|) with the expansion of 
the shape function operator, Eq. I|29|l . They are given by 



Wi(n-p,n-p = p + ) = S(u 2 ,m b ) 



(35) 



a s C F 

47T 



\n(p+ / p 2 )9{p+) 



l'2 



41n^-7 



l'.2 



-3 In 



9(p+) 



2 In 



2 n-p 

M2 



n ■ p 

M2 



,(/) 



where S{p,2,mb), and b± are the same as in section 

EH 
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-(k ■ v) = i{k ■ v)SZ, 

on c 

so that the wavefunction graph cancels Eq. (|34|l . when 
multiplied by the tree-level matrix element 

(k\O m \k) = (n-k) rn , 

I a — 5Z(n ■ k) rn = 0. The net result is that the opera- 
tor b v (in ■ D) m b v has no anomalous dimension from the 
graphs in Fig. at leading order in I/toj,. 

While this gives the correct result for the running of 
the local operators O m for m < 3, it does not repro- 
duce the known result for m > 3. For the latter case, 
however, there are additional graphs such as those in 
Fig. JHJ, which are divergent for m > 3 and renormal- 
ize the local operator. These box graphs connect the 
heavy quark to a light quark, and mix into four fermion 
operators. Such graphs are absent for deep inelastic scat- 
tering, where there is a twist expansion. Fermion fields 



VII. CONSTRAINTS ON THE ANOMALOUS 
DIMENSIONS 

There is a relation between the anomalous dimensions 
in the effective theory similar to the one obtained for 
deep inelastic scattering |19|. It is simplest to write the 
relations in terms of moments. The theory below /Xi has 
SCET currents, and depends on the label momentum 
fi ■ p. Anomalous dimensions can depend on n • p, but 
not on any infrared scales in this theory such as n ■ p, 
and so cannot depend on the moment N. In the theory 
below p2 j the effective theory is written in terms of shape 
function operators that depend on the scale p + = mb(l — 



2 The Feynman integrals for Fig. ||J for deep inelastic scattering 
involve fyfi = 0, whereas for B decay involve fqb ^ 0. 
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x) = rrib/N. The scale p~ has been integrated out, so 
the anomalous dimension is only a function of m^/N. 
Finally, in the theory below /13, the scale p + has also 
been integrated out, and the theory only knows about 
N, since the local operator depends on the moment. The 
anomalous dimensions 7* in the three theories below fii 
can be written as 



271 (aO = / 



JL 

n ■ p 



72(a) = / ( — ,o:«(m) ) 1 

73 (M) = f{N,a s (fi)). 

We have used twice the anomalous dimension 71 of the 
SCET current below fi\, since the time-ordered product 
involves two currents. 

The matching conditions at fi2 depend only on the in- 
variant mass p 2 of the hadronic final state which is in- 
tegrated out, so the matching coefficients are functions 
C2(p 2 /m 2 ) = C 2 (m b n ■ p/(Nfi 2 )). The matching con- 
ditions at /13 depend only on the scale p + which is in- 
tegrated out, so the matching coefficients are functions 
Cz(jp + /n) = C 3 (rrib / {N fi)) . Finally, using the relation 
that the derivative of the matching coefficients with re- 
spect to fj, is equal to the difference of the anomalous 
dimensions on cither side, one finds that the anomalous 
dimensions: (1) must be linear in ln^i, and (2) must have 
the form 

271(0) = A(a.(n))to.Jt-+B 1 (a a (n)), 
n ■ p 



72(0) = -A(a s (n))ln— +B 2 (a s (n)) 
73(0) = 0hiN + B 3 (a s (fi)), 



(36) 



which restricts the coefficients of the ln/x terms. The 
results Eqs. (|10I55|) satisfy these constraints, with A — 
—2oi s Cf I Ti"- The third relation shows that the anomalous 
dimension of b v (in ■ D) m b v should not have a lnm term 
as m — > 00. 

The anomalous dimensions 71,2 are used for fJ>i > [!> 
/i2 and fi2 > fJ- > 03, respectively. The ln/i term in 71 
is not a large logarithm at the upper endpoint fit, and 
the In fi term in 73 is not a large logarithm at the lower 
endpoint /13. At these limits, the effective theories match 
with full QCD above fit, which has no ln/i anomalous 
dimension, and the theory with local operators below /i 3 , 
which has no In N anomalous dimension. 



A. Relation to the factorization form of the 
anomalous dimension 

Korchemsky and Sterman j2^| have derived an alter- 
native form for the running of the B — » X s j moments 



between the scales nib and rrib/N. They have a renor- 
malization group evolution of the form 



c (7) = c ^ eA > 



where 



A 



1 dy 

t/N v 



T(a(m b y)) + ■y(a(m by /y)) 



-2T c (a(fi)) 



(37) 



Changing the order of integration gives 



A = 



m/VN A* 

m b /N 



In ^2T c {a(fi)) + T(a(fi)) + 2 7 (a(/x)) 
A* 



dfi 
fi 



In ^-2T c (a(fi)) + T(a(fi)) 



and so is equivalent to integrating the anomalous dimen- 
sion 271 between nib and nn/VN and the anomalous 
dimensions 72 between rrib/vN and nib/N where 

2 7l = -21n-^-r c (a(M))+rK/i))+27K/i)), 
m b 

72 = 21n^r c (Q( M )) + r(a( A( )), 
m b 

and so satisfies the constraint on the coefficients of the 
ln/x terms in Eq. (|36[) . with T c = u s Cf /n- Alternatively, 
the integrations from fi\ — > /i 2 and fii — > fi 3 can be com- 
bined into the single integral Eq. I|37|) since the constraint 
Eq. I|36fl is satisfied. 

The double integral form Eq. (|37|) has a Landau pole 
singularity for large moments. As for deep inelastic scat- 
tering, this Landau pole singularity is resolved by the 
effective theory computation [13. For large enough mo- 
ments, the scales /12 = rrib/VN and /13 = nib/N get 
smaller than /14 = Aqcd, and one needs to stop the 
renormalization group evolution at fi^. For B decays the 
condition [13 > /i4, which is required for Eq. H37J) to be 
valid is only true for the first two moments. 



VIII. APPLICATIONS 
A. Moments of the shape function 

It is generally believed that moments of the structure 
function f(k+) are related to local operators in HQET. 
This relation is used to construct models of the structure 
function using input from matrix elements of local op- 
erators extracted from inclusive B — > X C £D decays or to 
constrain fits to the shape function measured from data. 
We have already seen that this relation cannot be cor- 
rect, since the shape function and local operators have 
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different anomalous dimensions. In this section we show 
in more detail why this relation is incorrect. 

Moments of the bare operator oi°^ (k+) are given by 
bare local operators, 

/oo 
dfc + (-k + ) n Oi°\k+) 
-OO 

= b v {in-D) n b Vl 

but the corresponding relation does not hold for renor- 
malized operators. The infinite moments, over — oo < 
r + < oo of the renormalized operator O v (p + ) are singu- 
lar. One way to see this is to take the infinite moments 
of the renormalization group evolution Eq. I)23|l . The in- 
tegral over j(r + ,£ + ) does not converge at r + = oo, so 
the infinite moments require additional regularization. 

To gain more insight into the running of these moments 
without encountering the divergence, one can also study 
the evolution of half- infinite moments, defined by 



M+(0) = 
M N (0) = 



dr^ 



.N—l 



.N-l 



O r H 



O 



Consider the half infinite moment M N (O). We find 



dr+ (r+) 

i 

dt+O (t 



N-l 



d£ + -f(r+,£ + )0(e i 



dr H 



.N-l 



(O) has the diver gence at r + — > oo noted earlier. 
For moments of the shape function only the half infi- 
nite moment (O) is relevant, since f(r + ) vanishes for 
r + < —A. Thus, moments of the shape function are not 
defined without an additional renormalization prescrip- 
tion and are not directly related to matrix elements of 
local operators once radiative corrections are included. 
Models for the shape function should therefore not be 
constrained to have width a 2 = —Xi/(3m 2 ) for exam- 
ple. The parameters of a given model have to be kept 
arbitrary and determined directly from a fit to the data. 

Since the shape function is not equivalent to local oper- 
ators at order a s , the moments of the shape function are 
not given by the matrix elements of the local operators. 
Instead, one has to compute the matching correction onto 
the local operators, as in Eq. (|29|l . In particular the first 
moment of the shape function (because of the definition 
of the moments with powers of N — 1, the first moment 
is the normalization of the shape function) is not unity. 
The matrix element of Oq = b v b v is unity to all orders 
in a s and leading order in 1/rrib, by heavy quark sym- 
metry, so the first moment of the shape function is given 
by the first moment of Vq, the matching coefficient onto 
the lowest order operator Oq- £>o has a non-zero first 
moment, as can be seen from Eq. (I61f) . 

Note that low moments of the photon energy spec- 
trum in B — > X s ^f are still given by the well known 
matrix elements of local operators Ai, p\ etc. This is 
because for these moments one integrates over the entire 
region of phase space and the traditional OPE is there- 
fore applicable. High moments of the photon spectrum 
no longer match onto local operators, and are discussed 
in Sec. lVTTTCI 



Since j(r + ,£ + ) = unless £ + < r + , and r + < in the 
integration region, we can restrict the £ + integration to 
-oo < £+ < 0. Let r+ = z£+: 



»±M- N (O) 

OLsCj, 



d£ + (t 

DO 

dzz"- 1 



.N-l 



5(1 



a s C F 



l + 21n 



1-2 

d£+ 



,N-1 



^i 



fi 3 



(38) 



This shows that the half-infinite moments (O) 
are not multiplicatively renormalized, because of the 
ln(— £ + / n) term in Eq. (|38|l . The half-infinite moment 



B. Expressions for differential decay rates 

Using the results of this paper we can obtain expres- 
sions for differential decay rates to order a s . In this sec- 
tion we will give results for the photon energy spectrum 
dT/dx-y in the limit — > 1 for the decay B — > X s j and 
the double differential decay rate dT / (dzdxe) in the limit 
xg, — * 1 for the decay B — ■> X u lv. For both decays we 
give the results in the region rrifc m&(l — ^> Aqcd 
where an expansion in local operators is valid, and in the 
shape function region mb(l — ~ Aqcd- 



B 



X s j for m b > m b (l — a 7 ) > Aqcd 



In this region of phase space, the differential decay rate 
can be expanded around the limit x 7 = 1 , but we are not 
in the shape function region and the non-perturbative 
effects arc still obtained by matrix elements of local op- 
erators. The photon energy spectrum is given by 



1 dr s 



(x 7 ) = Wj (n ■ p,p + = TOfc(l — x 7 )) (39) 
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where W 1 was defined in Eq. (|15|) . Using the expression 
of W as a convolution of a Wilson coefficient C 7 with the 
bilocal operator O v , (|2"5f together with the expansion of 
the bilocal operator in terms of local operators given in 
Pty wc find 



= J dr+C 7 (ra fc (l - Xry) - r+, fi 2 ) 



^P m (r+,/i 2 ) 



(B\O m \B) 
2m b 



(40) 



In writing this expression we have ignored the running 
between the scales /x 2 and [13 and are therefore per- 
forming the matching onto the bilocal operator and the 
matching onto the local operator at the same scale. This 
allows one to give an expression for the decay spec- 
trum dr/dx 7 directly, rather than for its moments. The 
complete expression for the moments is given in Sec- 
tion Ennii 

Keeping only the lowest dimensional operator 
(B\Oq\B) — 2m b and using the explicit expressions for 
the Wilson coefficients C 7 and T>o given in Eqs. (|26|l . 
(123, JM| and Q]J we find 



ol s Cf 

47T 



(42) 



up to terms E(x 7 ) which have vanishing moments as 
N — > oo. The scale at which a s is evaluated is not deter- 
mined by this fixed order result. 

A comparison of Eq. i|42|) with the complete order a s 
expression 8] shows that the terms we have computed 
agree, and the remaining terms are 

E(x) = 7 + .x-2a; 2 -2(l + x)ln(l-j;), 

with moment 

iV 2 - N - 4 



M N [E(x)\ 



AN + 2 7 
N(N + 1) N+ N 



(N + l) 2 (N + 2)' 



which vanishes as N 
number 



oo. Here ffjv is the harmonic 



JV 



i dr 



S(fi2, m b )l 6(1 - x 7 ) 



(41) 



47T 



ln(l — x 7 ) 



7/ + 



l n 2 4 + 51n4-5-^-)«5(l-^ 



HIT. 



>»T. 



1/ 



where . . . denotes perturbative terms of order a 2 as well 
as terms of order a s which are not singular 
Furthermore, we have not given explicitly the power sup- 
pressed terms proportional to matrix elements of higher 
dimensional operators, which can be included by retain- 
ing the operators O m> \. The explicit logarithms of 
ji-i I nib in Eq. I|41|) are canceled by the [12 dependence of 
S(/i2, TTib). To fixed order in a s our results agree with the 
known perturbative corrections (see the next subsection), 
but we improve this result by including the Sudakov log- 
arithms originating from running between the scales m b 
and fi2 ~ 1-5 GeV. 



2. The B — > Xs'y spectrum to order a s ignoring 
renormalization group evolution 

The photon spectrum in B — * X s ^ decay to order a s 
ignoring renormalization group evolution is given by us- 
ing Eq. I|4ip. and expanding S(fi2,mt,) to order a s : 



1 dr 



ln(l — x 1 



1 — X* 



4tt 



2 1 



6(1 — x 7 



B 



X s -y for m b (l - x 7 ) ~ Aqcd 



In the shape function region, the expansion of the bilo- 
cal operator in terms of local operators can not be per- 
formed, and wc arc left with 



W-y 



rfr+C 7 (p+-fc + ,/i 2 )/(fc+,/i 2 ) (43) 



where p + = m b (l — XrS). Using the expression for C 7 
given in Eq. (|26|) we find 

1 dr s 



Tq Xj 



where 
G 8 (w,n) 



m b S((i2,mb) 
| a s C F 

47T 



V 2 ^ 



f(p + ,^) (44) 
dk+G s (p+ -k+, l i 2 )f(k+, l i 2 ) , 



2 
+4 



-5- 



7n 2 



S(w) 



In ^ 



- 21n^ 



311- 

U) 



(45) 



Note that this expression does not agree with previ- 
ous prescriptions to obtain the photon energy spectrum 
in the shape function region [l4|, where the perturba- 
tive spectrum is evaluated with a shifted b quark mass 
mb — * mb + k + and k + is then convoluted with the shape 
function f(k + ). The reason for this disagreement is dis- 
cussed in more detail in Section EnHH 
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4- B — » /or m 6 > m(,(l - x e ) > Aqci 



where p + = m b (l — xp), 



Any differential decay rates can be obtained from the 
general expression of the triple differential rate given in 
Eq. ©. As an illustration, we present results for the 
double differential decay rate dT u /(dzdxg) in the limit 
xi — * 1. Using p 2 <g; 1, we find for this double differential 
rate 



1 dr u 



= 12^(1 -z) f 
Jo 



dp+ 



2Wi + m b zW%] . 



(46) 



As for the decay B — > X s j, one can use the expression for 
as a convolution of a Wilson coefficient cj with 
the bilocal operator O v , Eq. 125( 1. together with the ex- 
pansion of the bilocal operator in terms of local operators 
given in Eq. (|2"5|) . Using the results given in Section IVEl 
for the Wilson coefficients Cf together with Appendix lAl 
to perform the integration over the p. distributions in the 
Cf we find 



1 dP 1 
Y% dzdxt 



12z(l- z)S{p 2 ,m b ){ 1 + 



a s C[, 

47T 



with 

H(z,xi) 



2 2 

In 2 t 2 u ln z _ 5 x ln jh_ _ 2H(z, xi) 

m b m b 



JhYl -.,-,)- [2\nz- 7 -)\mi -- .,-,) 



-2Li 2 (l-z) 



(3z- l)lnz 
2(1 -z) 



27^ 

3 



',(47) 



In 2 z 



(48) 



As for the B — * X s j case (section I VIII B 111 we have ne- 
glected the running below the scale p2 and the depen- 
dence on H2 in Eq. I|47ll cancels. Expanding S f (/i2,m; ) ) to 
first order in a s we reproduce the result obtained previ- 
ously obtained by De Fazio and Neubert |14| . 



B 



X u lv for m b (l - x e ) ~ Aqcd 



As for the decay B — » X s j the expansion of the bilo- 
cal operator can not be performed in this region of phase 
space and differential decay rates in this region of phase 
space are given in terms of a convolution of a perturba- 
tively calculable coefficient and the shape function. Us- 
ing the results of Section IV El for the Wf in the shape 
function region we find 



1 dr u 

Tq dzdxi 



12z(l - z)S{p2,mb) 



F(p+,p 2 ) (49) 



47T 



f 'dk + G u (p+ -fc+, M 2)/(fc + ,M2) 
Jo 



F(p+,p)= / dfc+/(fc + ,/i), 
Jo 



(50) 



and 



G u {uj 7 fi) — 21n 31n 4 In . 



-4Li 2 (l-z)-21n^^-5-^-. (51) 
1 — z 6 



Again, our results disagree with the results presented in 
Ref. as we will discuss next. 



6. Comparison with previous results in the literature 

De Fazio and Neubert include shape function effects by 
convoluting the perturbative decay spectrum (i.e. Wf 3 **) 
with a shape function. This procedure is not correct be- 
cause there are large perturbative corrections included 
in the definition of the shape function. Consider, for ex- 
ample, the most singular terms in Wf art , the lnp + /p + 

terms. These arise from (p + , p) + /p2t(P + : A 1 )- The 
term gives a coefficient of 4, and the f^J^ term gives 
a coefficient of —8, for a net contribution of —4. The 
order a s computation of Ref. gives the net result of 
—4, without reference to the scale at which it arises. Our 
effective theory computation shows that it is generated 
at two different scales, 4 is generated from the match- 
ing condition at P2 onto the shape function operator and 
—8 from the matrix element of the shape function at 113. 
De Fazio and Neubert convolute the shape function with 
the order a s spectrum, i.e. with the —4. Instead, one 

should convolute it with the order a s coefficients C\ ' 
generated by the time-ordered product of currents at the 
scale H2, which has a coefficient +4. The remaining —8 
is included in the definition of the shape function, since 
it is produced by taking the matrix element of the shape 
function operator in the B meson. This somewhat sur- 
prising conclusion follows from the fact that the the shape 
function operator has a non-zero one-loop matrix element 
between on-shell b quark states, as shown in Eq. I|31|l . 
In deep inelastic scattering, where the parton distribu- 
tion operator has zero matrix element between on-shell 
quark states, the cross-section can be written as the con- 
volution of the order a s perturbative cross-section with 
a non-perturbative parton distribution function. 

Below the scale 113, the entire order a s coefficient has 
been removed from from the operators, and included 
in the coefficient functions. The theory below /13 is 
an expansion in local operators, which is valid when 
p + Aqcd, and so is no longer the shape function re- 
gion. 
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C. Renormalization group improved photon energy 
moments 



Moments of the photon energy spectrum as defined in 
[ll2^| are given by (using Eqs. @ED13) 



m B /m b i j-p 

1 ax j 

m B /m b pm b (l—x 7 ) 

dx-yx^' 1 / dr H 



C 7 (rrib(l — £ 7 ) — r + , /i) 

rm B /m. b 



m b 



(B\O v (r+)\B) 
2m b 

m B /m b 

dy 



C 7 (m b (y - 



(B|0 B (m 6 (l-i/))|B) 



2m? 



where we have made the change of variables r + = m b (l — 
y). Note that the physical limit on the photon energy is 
ms/2, making the upper limit on x 7 to be ms/mb- The 
limits of integration on r + arise because the Wilson coef- 
ficient vanishes for r + > m b (l — x 7 ) and f(r + ) vanishes 
for r + < —A. Changing the order of integration, and 
making the change of variables x 7 = yz gives 



Ml = 



m B /m b 



dy y 



N _ 1 {B\O v {m b (l~y))\B) 



2m b 



< I dz z 77i b yC~f {m b y(l - z), (i) . 
lo 

The coefficient C 1 (p + ) has terms of the form S(p + ), 
[l/p + ] ll and [\np + /p + ]^. Using the results of Ap- 
pendix [S] Eqs. 1A2IA3H one can convert to the form 

C 7 (m b y(l — = C 7 (l — z, fj,) + hiy terms, 

where C 7 (l — z) contains 5(1 — z), [1/(1 — z)]+ and [ln(l — 
z) / (1 — z)] + terms. In the endpoint region y — > 1, lay = 
In [1 + (1 — y)] — > 0, so the lny terms can be dropped. 
Thus Mn is a product of moments, 



Ml 



M N (f)M N (C 7 ), 



(52) 



where the moments of the shape function and the C 7 are 
defined by 



M N (f) = 
M N (C 7 ) = 



m B /m b 



dy y 



N 



_ 1 (B\O v (m b (l-y))\B) 



dz z Ar - 1 C 7 (l-z, M ). 



2ms 



The analog of Eq. (|52"|> holds for any , with C 7 re- 
placed by C { /\ 



M 



N (Wf } ) = M N (f) M N (cj /} ) (53) 



where 
M N (wty 



l/m b 



dx x N 1 W(n ■ p,p + = m b (l — x 
The moments of the shape function operator 

rm B /m b 

N (fi) = / dvv tf - 1 O v (m b (l-v)), 
Jo 

where 



M N (f) = 



(B\O n (vl)\B) 
2ms 



satisfy the renormalization group equation for its matrix 
element 



s~i pm B /m b poo 

Of / ' N _ 1 
dyy 



dx (O v (m b (l - x))) 



xi5(x-y) 
+2 



((x-y)m b >0 +ln j §{x _ y) 

x-y (i 



using the anomalous dimension Eq. I|24|) with r + = 
m b (\ — y) and t + = m b (l — x). Letting y = xz gives 

Mj" (Oat (^)) 
d/x 

n (In r^Blm b f -l 

■ r 1 dxx N - 1 (O v (m b (l-x))) dzz"- 1 



xi 1 - 2 In 



m b x 



5(1 -z) + 2 



1 - z 



Since we are looking at large moments, x — * 1 and In a; — ► 
0, so that 

M A( 0jv ( M )) = _ 7jv(// ) , (54) 

d// 

where the anomalous dimension 7jv(a*) is, at one loop, 



7jv : 

which for large iV is 

7at 



l_21n-^-) -2H N -n >. 
m b 



^— - 1 - 2 In — 

7T V mb 



(55) 



This agrees with the results of The complete NLO 
anomalous dimension is given by including the two-loop 
In N part of the anomalous dimension, which is given by 
the relation Eq. (|36|) . 



In = 45 



(27T) 



TOfc 



(56) 
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where B is given in Eq. 1111) . 

The solution to this renormalization group equation is 
given by 



(0 N {V3)) = S' / 1 (M3,M2)(Ojv(M2)> 

where the A dependent scale factor is 

r fo(z) 



exp 



+ r f i(z) 



(57) 



(58) 



with 



r fo(z) 
r fi( z ) 



8wC f 



PI 




2C F p 1 


Pi 




2C F 


2 


Po 




4C F B 



1 + In z 

z 



1 — z -\-\nz — — In z 



In ; 



where 



z — 1 — In z 



(59) 



(60) 



Combining these results we obtain for the renormal- 



ization group improved moments of W i 



(/) 



M 



N 



(/) 



M 



N 



cfW 



Mjv[/(/x 2 )] Mjy [cfV 2 ) 
5( M2 ,Mi) a(%w[CK/ia)] 



q s (m b )C F &(/) 

47T ' 



Mjv [C(w,/i 2 )J = H Mtv 



4tt 



W,M2, 



Mjv[/(/i 2 )] = 5/(M3,M2)A^[/(M3)]. 

At the partonic level, 

AM/( M3 )] = l + ^^^M^M]. 



-in 



The moments of the pcrturbative coefficients are, using 
Eq. JIU, 



N-l 



M 



N 



CWK M2 ) = 4 ^^-(41n^^-:!)//A , 



/'2 



21n2 m^_ 31n m^ + 7 _ 7r2 



Ma 

7V-1 



M2 



M W [/W( M3 )] = -8^^-4(ln^-l)^_ 1 



ln 2 4 + 21n4-^ 



m 2 6 



(61) 



These expressions can be used to obtain the renormal- 
ization group improved moments for the B — > A s 7 spec- 
trum. Using = 1, 



M N 



i dr 



W, 



x<M + 



47T 



a s (m b )bW + a s (fi 2 )M N [c (1) (w,/x 2 ) 



a s (mb)Cf 

47T 



Mat [JS] 



(62) 



The computation to leading order in the SCET expan- 
sion parameter A does not include the terms E(x) with 
vanishing moments as A — » oo. We have included these 
terms so that we can compare the renormalization group 
improved moments with the fixed order result even for 
small moments. Since E(x) is generated in the tradi- 
tional operator product expansion at the scale nib, we 
have included these terms with a s (rrib). We have not 
included the running of local operators below /i 3 . As dis- 
cussed below, the regime 1x3 > [i^ exists only for the first 
few moments, and for these, the local operators have zero 
anomalous dimension. 

We can now compare the order a s moments with the 
renormalization group improved moments. For simplic- 
ity, we will compare the two parton model results, i.e. we 
will use Eq. H61fl for the moments of the shape function. 
Non-perturbative effects can be included in Mm{J) if so 
desired. We choose /ii = mj and /14 = 1 GeV, which 
are fixed. /12 = to^/VA and ^3 = nib/N vary with the 
moment. For A > 3, /i 3 < fi£. In this case, the regime 
of local operators does not exist. One runs the shape 
function down to the scale ^4, and computes its matrix 
elements in the B meson state at /Z4, so one can effec- 
tively set /13 = /14 in Eq. i|62|) . For even larger moments, 
A > 14, fi 2 < /M- In this case, the shape function regime 
does not exist. One runs the SCET currents down to 
Hi and then takes matrix elements of the time ordered 
product in the B meson state. In this case, one can ef- 
fectively set H2 = = ^4 in Eq. i|62|l . The order a s 
and renormalization group improved moments are show 
in Fig. \7\ The renormalization group resummation pro- 
duces a significant correction to the moments, and cannot 
be neglected. 

The large moments in the perturbation theory calcu- 
lation shown in Fig.[7]become negative. This a reflection 
of the breakdown of perturbation theory. In obtaining 
Fig. [3 we have stopped the renormalization group evo- 
lution at the scale /xj, but still computed the B matrix 
element using perturbation theory. The matrix elements 
in Eq. H61|l contain In A terms, which are an indication 
of the sensitivity to the lower scales nib/VN and m^/A, 
and cause a breakdown of perturbation theory for large 
N. The correct computation for large moments involves 
taking the non-perturbative matrix elements in the B 
meson state. In the shape function region (/X3 > /X4 but 
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FIG. 7: The first twenty moments Mn of the B — ■> X a y pho- 
ton spectrum (l/ro)(dr/da: 7 ) — renormalization group im- 
proved (solid black), order a 3 (dotted blue), and tree level 
(dashed red). The kink in the blue curve at N ~ 2.75 is the 
point where ^13 = {14, and the kink at N w 13.5 is where 

M2 < ^4)7 one sets 113 — [i^ in Eq. (|62|l and uses the 
non-perturbative results for the shape function moments 
A/jv[/(M3)]- For /i2 < /i4, one sets P2 = A*4; drops the per- 
turbative corrections generated at the scale and 
uses for Mjv[/(a*3)] the moments of the non-perturbative 
matrix element of the time-ordered product of currents. 



1. The decay spectrum vs. its moments 

The effective theory computation is most naturally 
done in moment space. The intermediate scales \i2 and 
/X3 are given in terms of the moment N as mt/VTV and 
m,b I N. The matching contributions and running depend 
on /i2 and ^3, which change with N. Furthermore, for 
large N, the scales ^2,3 can become smaller than Aqcd, 
and the effective theory computation has to be modified, 
as discussed above. 

We have not found a simple way to give a theoretically 
correct expression for the spectrum in x-space which in- 
cludes the renormalization group evolution, and at the 
same time properly takes into account effects such as fj.2,3 
becoming smaller than Aqcd- The simple double inte- 
gral form Eq. I|37(l is not valid because it has a Landau 
pole singularity as x — ► 1. 

IX. CONCLUSIONS 

We have calculated the hadronic tensor for left handed 
vector and tensor currents, relevant for the decays B — > 



X u lv and B — * X s 7 in a region of phase space where 
the final hadronic system satisfies p + <C p~ . From this 
hadronic tensor any differential decay rate can be ob- 
tained. This calculation has been performed using the 
soft collinear effective theory, by performing the match- 
ing in multiple steps and summing all the logarithms aris- 
ing from the different scales. We considered two possible 
regions, where p + ~ Aqcd and where p + 3> Aqcd- 

In the former region, the hadronic tensor is given as 
a convolution of calculable Wilson coefficients and the 
lightcone distribution function of the B meson. We cal- 
culated these Wilson coefficients to order a s including the 
leading and subleading Sudakov logarithms. The result- 
ing expressions, when expanded to order a s (i.e. with- 
out renormalization group improvement), agree with the 
order a s results previously given by de Fazio and Neu- 
bert J14|. However, we do not agree with the assumption 
in Ref. |14| that the structure function of the B meson 
does not contain any perturbative effects, so that shape 
function effects can be included by convoluting the or- 
der a s decay distributions with the shape function. We 
have seen that very large perturbative corrections are in- 
cluded in the shape function matrix element Eq. Li 21) at 
the scale ^3. We have also shown that the usual assump- 
tion that moments of the lightcone distribution function 
are given by matrix elements of universal local operators 
(such as the kinetic or Darwin operators) is incorrect. 
This relation only holds at the classical level, but fails 
once quantum corrections are taken into account. 

In the region p + 3> Aqcd the structure function can 
be expanded in terms of local operators and the the 
non-perturbative physics of the hadronic tenor are pa- 
rameterized by matrix elements of these universal local 
operators. The results presented here contain the most 
singular contributions in the endpoint region, which are 
enhanced by powers of m,b/p + relative to terms we have 
dropped. The matching coefficients for these enhanced 
non-perturbative contributions are again given to order 
a s . 
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APPENDIX A: THE fi DISTRIBUTION 

It is convenient to use a modified version of the usual 
+- distribution to regulate singular functions of p + on the 
infinite interval p + S [0, 00]. The /i-distribution is defined 
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by 



e(p+) 



ln(p+//i) 6(p+) 



: lim 



= lim 

5^0 



^W)ln^ 
p+ /i 

ln(p+/ M ) 6(p+ > 



2 £ 



--5(p+)ln 



and is closely related to the ^-distribution in Ref. [U 
The /i-distribution satisfies the relations 



-A 



dp + f(p+) 



6(p+) 



dp+ /(p + )-/(Q) 



/(0)ln- 



D 



dp+ f(p + ) 



m(p+//z) %+)" 



(Al) 



for A > and £> > 0. The /^-distribution is defined 
so that integrals such as Eq. i|Al[) involve logarithms of 
dimensionless quantities. 

The /i-distribution satisfies the scaling property 



9(Xp+) 
Xp+ 

H\p+/p) e(x P +) 



Xp^ 



9(Xp + )] In A 
m(p+//z) 9(p+) 



In A 



0(Xp+) 



Ap^ 



In 2 A 



+ (A2) 



One can also consider the /i-distribution restricted to 
a finite interval p + € [0, A]. In this case, one can convert 
to the usual -{—distributions defined by 



,A 

/ dp+/(p-» 
Jo 





-/ 




+ Jo 



dp+ /o+) 



ln(p+/ M ) 9(p+) 



B 



dp 1 



/(p+)-/(0), p+ 



/o /' 
The conversion is 

r%+) 
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hi: 
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ln(p+/ M ) 9(p+) 



\6(p+ 
[ P + 


-1 + S( P + 

- + 


"ln(;U 


7m) o(p + )] 


"ln(^ 


p+ I 

7A) 0(p+)"| 


L p + J 



A 



InA 



%+) 



i<5b+)ln 2 ^ 



^( P + )ln 2 ^, 
2 // 



(A3) 



which depends on the size of the interval A. We will need 
this conversion formula with A = m& andp + = m^l — z), 



6(p+) 



1 



Hp+/fi) 6(p+) 



1-z 
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Using this, Eqs. (|17I32|I become 
CW(g+, M2 ) = 4 
2 In 



"ln(l - z)~ 




[ (1-z) _ 





V-2 



1 - Z 



- 3 In ■ 



/« (fc+, M3 ) = -8 



ln(l - z) 



— ^ + 7-7T 2 U(l-Z), 

-l|ln4-l, , 

mt I 1 — z 



+ (-ln 2 4 + 21n4-^)^ | --» 



(A4) 
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